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Abstract. The spin-down of Be stars due to angular mo- 
mentum transport from star to disc has been considered. 
This has been prompted by empirical studies of observed 
optical and IR line profile studies indicating that the disc 
is rotating in a Keplerian fashion. It is found that sub- 
stantial spin-down may occur, especially for late B stars 
throughout their main-sequence lives for the "strongest" 
discs (most dense ~ 10~^^g cm~^ with high radial velocity 
^ 1 km s^^ at their inner edge and with large opening an- 
gle ~ 15°). This is in conflict with studies of rotational ve- 
locity distributions for different luminosity classes, which 
show no significant evolution. The implications of this for 
Be star discs are considered. 

Key words: stars: emission-line. Be - stars: rotation - 
stars: evolution - circumstellar matter 



1. Introduction 

Be stars are now widely accepted to have two distinct 
regions of circumstellar matter : a diffuse polar stellar 
wind and a dense equatorial "disc" (Dachs 1987, Slette- 
bak 1988). One of the major objectives of Be star research 
is to develop a theory which describes both of these com- 
ponents. The fast diffuse polar wind is well described by 
standard line-driven wind theory (Castor, Abbott & Klein 
1975, Friend & Abbott 1986, Kudritzki et al. 1989). How- 
ever, it has been much more difficult to describe the equa- 
torial disc. Empirical models of the disc structure have 
been presented by e.g. Marlborough (1969), Waters (1986) 
and Hanuschik (1996), whilst theoretically driven models 
have been developed by Poe & Friend (1986), Chen & 
Marlborough (1992), Bjorkman & Cassinelli (1993), Will- 
son (1986), Ando (1986) and Lee, Saio & Osaki (1991). 
These involve phenomena such as magnetic winds, latitu- 
dinal variation of driving lines, wind compression, stellar 
pulsation and viscous excretion. The most promising disc 
theory for several years, Bjorkman & Cassinelli's (1993) 



wind compressed disc model, has been shown to be inca- 
pable of reproducing observed discs by Owocki, Cranmer 
& Gayley (1996) and Porter (1997) via different routes. 

One feature of empirical studies of line profiles in Be 
star discs is that they imply a rotationally supported disc, 
and that the rotation falls off in an Keplerian fashion (e.g. 
Dachs et al. 1986, Hanuschik 1989, 1996). The half-line 
width is typically larger than usini (Hanuschik 1996). 

Also, the current model for V/R variations (described 
in e.g. Dachs 1987) assumes that m — 1 perturbations 
arise in a Keplerian disc (Papaloizou et al. 1992). Excre- 
tion disc models proposed by Lee, Saio & Osaki (1991) 
provide naturally this sort of rotationally supported disc 
- the disc material is rotating at its Keplerian speed, and 
drifts outward due to the effect of viscosity. These models 
require that angular momentum is supplied at the inner 
boundary of the disc. Given a prescription for the viscosity 
(e.g. 'a' from Shakura & Sunyaev 1973), then the surface 
disc density, and disc scale height may be integrated from 
the equations conserving angular momentum and mass 
(e.g. see Pringle 1981). A similar model has been used by 
Pringle (1991) applied to the cessation of accretion ("de- 
cretion") by a young stellar object. 

In this investigation, assuming that the disc is indeed 
rotationally supported, the spin-down of the central star is 
calculated. In §2 the evolution of a star's rotational veloc- 
ity is derived given that it is supplying angular momentum 
to the disc. Estimates of spin-down times are presented in 
§3 across the B star range of stellar parameters. This is 
discussed in §4 and conclusions are presented in §5. 

2. Angular momentum transfer & spin down 

It is here assumed that the disc material is moving round 
the star in approximately a Keplerian fashion. The star, 
however is not rotating at its Keplerian velocity, and so 
some angular momentum needs to be transferred to the 
disc. This may be achieved through non-radial pulsations 
(described by Osaki 1986) or through the action of mag- 
netic fields. However, the exact mechanism is not impor- 
tant for the discussion here. 
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The specific angular momentum of the gas at the star's 
surface U, rotating at velocity v^^* is 



7 - p _J2GM, 



(1) 



where R^ and Rp are the equatorial and polar radius of 
the star respectively (i?* > Rp), is the mass of the 
star, and the star is rotating at a fraction / of its critical 
(or break-up) velocity (G is the gravitational constant). 

This should be compared to the specific Keplerian an- 

1 

gular momentum 1^ = v^^^r = {GM.^ r) 2 . It is now as- 
sumed that the extra angular momentum is added to the 
disc at or close to the star's surface - this is a conservative 
estimate leading to a lower bound on the angular momen- 
tum deficit. If the disc has a half-opening angle 9, and 9 
is small, then the rate of angular momentum supplied to 
the disc is 



dL 



= (Ik - QAnR^'^ej- 



\^Rp 



(2) 



where J- is the mass flux through the disc in g cm^^s^^. 

The rate of angular momentum lost by the star is then 
simply dL^,/dt = —dL/dt. If the star rotates as a solid 
body at angular velocity and is regarded as a poly- 
trope of index 3/2, then the total angular momentum of 
the star is 



/2GRA 2 



V 3 



(3) 



where e = 0.2046 (this is a slight underestimate as the 
star's mean radius has been assumed to be Rp). This leads 
to the rate of angular momentum loss by the star: 



dL^, T, , 3 f 2GRp 

= eM*2 tL 

dt \ 3 J dt 



2 df 



(4) 



which may be combined with eq.2 to finally give 



df _ _ ( Ane^( TR^'^^(2,R 
di~~ 




(5) 



2.1. Disc parameters 

To calculate the spin-down, the parameters of the disc 
need to be fixed. Many of the models for Be star discs 
may be described via power law expressions of density, 
disc height and velocity. For example, both Lee et al.'s 
(1991) and Okazaki's (1997) excretion disc models may 
be represented by power law density distributions. 



The discs must also be normalized in some way such 

that they provide sufficient IR excess as the various disc 
models outlined above do not yield absolute values of den- 
sities, mass-loss rates etc. Here the density at the inner 
edge of the disc is used from Waters, Cote & Lamers' 
(1987) fitting to IR continuum excess. 

The mass-flux in Waters et al.'s model is ^ = 
7.7 X 10~^uo/0-ii where the density p-n is measured in 
10~^^g cm^'^, the velocity at the inner edge of the disc 
vo in km s~^and in g cm~'^s~^. Waters et al. calculate 
densities of w 0.05-2 (their table 4) across a range 

of B spectral sub-types. The velocity of the flow assumed 
by Waters et al. was 5 km s~^ (yielding a high mass flux) 
is almost certainly too high (Marlborough et al. 1997). 
Exactly what value of the initial velocity should be used 
depends upon which model is implemented, although vq 
has been calculated to be subsonic from both empirical 
and theoretical models. Poeckert et al. (1982), Marlbor- 
ough & Cowley (1974) & Marlborough et al. (1997) find 
initial velocities of ^ 1 km s^^ using empirical models, 
whilst for excretion disc models it can be shown that the 
radial velocity of an "alpha" accretion (or excretion) disc 
is Vr ^ acl/v^, where Cg is the sound speed (e.g. Pringle 
1981) which becomes 



Vo w 0.3ar4 



Ml 



km s 



(6) 



where T4 is the disc temperature in 10'*K, and and 
R^, are measured in Solar units. This leads to values for 
the velocity of vo = 0.2a-0.6a km s~^across the range of 
stellar parameters for stars (sec below) . The final unknown 
in eq.5 is the disc half-opening angle 9. This have been 
calculated statistically to be ~ 5° (Porter 1996) to ~ 
13° (Hanuschik 1996) for Be star samples, whilst Waters 
et al. (1987) use an opening angle of ^ = 15°. 

With the above mass flux normalized to a base den- 
sity, this should ensure that the disc should be able to 
produce enough free- free emission in the IR to account for 
the observed cxcesscis. However, this is still subject to the 
functional form of the density distribution and thickness 
of the disc model. The above normalization is therefore 
a necessary, but not sufficient procedure for all disc mod- 
els. The study by Marlborough et al. (1997) illustrates the 
effect of differing power law density distributions on the 
optical and IR line profiles. 

3. Results 

Combining the estimate for mass flux with eq.5 gives 



Rr, J \ \ "^Rp 



(7) 
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Fig. 1. The function f{t) given by eq.7 for different disc 
models and different stellar parameters. For each panel the 
solid, dashed and dash-dotted lines correspond to the f 7.5M0 , 
7.7Mo and 3AMq models. An initial value of / = 0.81 at t = 
is assumed. 



where Rp is measured in Solar units. This function may 

now bo integrated to obtain the evolution of the rotation 
of a Be star. Three model stars (with stellar parameters 
in Table 1) spanning the B spectral type are used. 

Eq.7 breaks into three separate parts: the stellar pa- 
rameters, the disc parameters, and a function of rota- 
tion. For a given star the models may be characterised 
by the parameter P = voOp-u. With 6 = 5-10°, 
p_j^^ = 0.05 - 2.0, and Vq = 0.1-1.0 km s-\ the disc 
parameter ranges between P = (0.5-350) x lO"*^. Three 
values are used; P = 5 x 10~^, 10~^, and 10~^ correspond- 
ing to "weak", "medium" and "strong" discs. 



Table 1. Stellar parameters for the stars considered from 
Schmidt-Kaler (1982). The radii quoted are the polar radii, 't 
is an approximation to the main sequence lifetime taken from 
the hydrogen burning times from Maeder & Meynet (1989) 



Spectral 


M. 






type 


(Me) 




(yr) 


BO 


17.5 


7.7 


9.8 X 10" 


B3 


7.7 


4.7 


3.9 X 10^ 


B9 


3.4 


2.8 


2.7 X 10® 



Figure 1 displays the spin-down due to each of these 
discs with the top, middle and bottom panels correspond- 
ing to the strong, medium and weak discs respectively. The 
initial value of / was set to 0.81, and the calculation for 
each star was stopped at a time equal to the main sequence 
lifetime. This was taken to be the time at which hydrogen 
burning ceases in the models of Maeder & Meynet (1989). 
These times are given in column 4 of table 1. 

The results displayed in fig.l clearly illustrate three 
things: (i) the amount of spin-down is strongly dependent 
on the disc parameter P, (ii) the amount of spin-down for 
a star is not strongly dependent on its stellar parameters, 
except through its lifetime, and (iii) for some combinations 
of disc and star, substantial spin-down will occur. 

For early B stars, there is only significant spin-down 
during the main sequence for strong discs. Both medium 
and weak discs cause little evolution in the rotation. How- 
ever, for late B stars unless P < 5 x lO"'' significant spin- 
down should occur. If there is a threshold rotational ve- 
locity below which no disc should form (previous studies 
have found this to be / ~ 0.2 — 0.3; Porter 1996), then 
late B stars may even spin-down so much that they cease 
to support a disc at all. 

4. Discussion 

It has been demonstrated above that there is a possibil- 
ity of some Be stars significantly spinning down. However, 
in observational studies, Fukuda (1982), Slettebak (1982) 
and more recently Zorec & Briot (1997) have foimd that 
across the B star range, there is no significant change in 
the rotational velocity distributions between different lu- 
minosity classes. This, contrasted with the above numer- 
ical results makes it very difficult to see how, (for late B 
stars), even medium discs may be generated, as these stars 
should have spun-down significantly during their main se- 
quence lives, or with strong discs they should have spun- 
down to small rotational velocities and may cease to sup- 
port a disc. As the disc structure has little or no depen- 
dance on the spectral type of the star (van Kerkwijk et 
al. 1995) then, at a first glance, it is perplexing that ob- 
servational studies indicate that the discs around Be stars 
rotate at Keplerian velocities. 

It is possible that all the actual Be star discs are what 
have been termed weak here. Some small fraction of stars 
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may indeed have strong discs and even though they may 
spin-down significantly during their lifetime they may not 
skew the rotational distributions enough to detect. The 
conclusion of this is that the disc population is dominated 
by weak discs. The least well known parameter in the disc 
is the initial radial velocity vq. Assuming maximal values 
of » 2 and 6 = 15° for the disc, then a star will 
not spin down significantly if the initial radial velocity is 
vo ^0.01kms-i(P ^ 5x10-3). 

The estimates of vo stem from optical and IR line mod- 
elling of discs although these only provide upper limits of 
~ 1 km s"^. Excretion disc models (e.g. Lee et al. 1991) 
do predict values of vq, although at the expense of intro- 
ducing the viscosity parameter a (cq.6). With a <^ 0.01, 
values of vq applicable to weak discs are obtained. Imposi- 
tion of no evolution in the rotational distribution therefore 
may be seen as a tighter constraint on vq (or on a). 

In the discussion so far, the disc has been assumed to 
exist around the star at all times in the main-sequence 
lifetime of the star. Be stars are known to undergo phases 
when the disc is apparently lost (sec e.g. Hirata 1995 and 
references therein). Depending on the fraction of their life- 
time spent without discs (or with very weak discs), Be 
stars could easily support strong discs and still not evolve 
significantly in rotational velocity. 

It is timely to discuss whether the disc material is gen- 
uinely being lost. Strong discs could exist around B stars if 
all, or a fraction of them are re-accreted at some phase and 
hence spinning-up the B star. Indeed, a similar process of 
viscous re-accretion of disc material has been suggested 
by Hanuschik et al. (1993) in their study of /z Cen. How- 
ever, at low radial velocities such as these it is impossible 
to detect whether material is inflowing or outflowing as 
the rotation dominates the kinematics in the disc (also 
see Hanuschik et al. 1993). This renders the line profiles 
generated in the disc impotent as radial velocity detectors 
in the presence of the Keplerian rotation. Consequently it 
seems impossible to obtain an explicit observational mea- 
surement of Vq. 



5. Conclusion 

Under the assumption that Be stars maintain discs 
throughout their main-sequence lives, spin-down is found 
to occur for certain parameters of the circumstellar disc 
if the disc rotates in a Keplerian fashion. This is in con- 
flict with observational studies of the distribution of rota- 
tional velocity for Be stars of different luminosity classes. 
The simplest resolution to this paradox is that the major- 
ity of Be star discs must occupy parameter space in which 
P = VQ0p^ii 5 X 10"^. For the range of opening angles 
from observational studies (e.g. Porter 1996), and densi- 
ties (e.g. Waters et al. 1987) this implies very low radial 
velocities of vq <^ 0.01 km s~^. 
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